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The rise of the cosmic ray positron fraction with energy, as first observed with high confidence by 
PAMELA, implies that a large fiux of high energy positrons has been recently (or is being currently) 
injected into the local volume of the Milky Way. With the new and much more precise measurement 
of the positron fraction recently provided by AMS, we revisit the question of the origin of these high 
energy positrons. We find that while some dark matter models (annihilating directly to electrons 
or muons) no longer appear to be capable of accommodating these data, other models in which 
~l-3 TeV dark matter particles annihilate to unstable intermediate states could still be responsible 
for the observed signal. Nearby pulsars also remain capable of explaining the observed positron 
fraction. Future measurements of the positron fraction by AMS (using a larger data set), combined 
with their anticipated measurements of various cosmic ray secondary-to-primary ratios, may enable 
us to further discriminate between these remaining scenarios. 

PACS numbers: 95.85.Ry, 95.35.+d, 97.60.Gb; FERMILAB-PUB-13-090-A 



I. INTRODUCTION 

Recently, the AMS collaboration published its mea- 
surement of the cosmic ray positron fraction over the 
range of 0.5 to 350 GeV [l|. Their findings confirm with 
unprecedented precision earlier measurements from the 
PAMELA [i| and Fermi [1| collaborations, which had 
each reported a clear rise in the positron fraction at en- 
ergies above ~10 GeV [4] (hints of such a rise were also 
present in data from HEAT [5,] and AMS- 01 [1]). At 
the time, PAMELA^s observation generated a great deal 
of interest and speculation as to the origin of the high 
energy positrons. Leading proposals put forth to explain 
this observation included dark matter (DM) particles an- 
nihilating or decaying in the Galactic Halo (7H21J. and 
nearby pulsars injecting high energy positrons into the 
interstellar medium |22H2a |. 

Despite providing valuable information, the measure- 
ments provided by PAMELA and Fermi were not suffi- 
cient to discriminate between DM and pulsar origins of 
the rising positron fraction. The much higher precision 
measurement of the positron fraction by AMS, however, 
brings new and important information to bear on this 
question. In this article, we make use of this new data 
and revisit both annihilating DM and pulsars as potential 
sources of the observed high energy cosmic ray positrons. 
We find that DM particles which annihilate directly to 
e^e^ or fJ.^ fJ-^ can no longer accommodate the observed 
positron fraction. However, DM particles with a mass 
of ~l-3 TeV annihilating to intermediate states which 
then decay to muons or charged pions could potentially 
provide a good fit. Pulsars also continue to represent a 
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potentially viable explanation for the observed positrons. 
The remainder of this article is structured as follows. 
In Sec. ini we discuss whether annihilating DM can ac- 
count for AMS^s measurement of the positron fraction, 
considering a variety of DM models and models of cos- 
mic ray propagation. In Sec. Illli we discuss whether 
pulsars can account for the observed data. In compar- 
ing these scenarios, we find that the existing data from 
AMS cannot yet definitively discriminate between the 
DM and pulsars origins of the observed positrons, al- 
though the range of models capable of accommodating 
the data is now significantly more constrained. With fu- 
ture data from AMS, providing not only measurements 
of the positron fraction but also of various cosmic ray 
secondary-to-primary ratios, we expect to be able to fur- 
ther narrow the range of models potentially responsible 
for the rising positron fraction. In Sec. lIVl we summarize 
our results and briefiy discuss our current understanding 
of the possible origins of the observed cosmic ray positron 
fraction. 



II. ANNIHILATING DARK MATTER 

If annihilating DM particles are to account for the ob- 
served rise in the cosmic ray positron fraction, they must 
be quite heavy, certainly no less than 350 GeV. Further- 
more, models which can also accommodate the smoothly 
varying and consistently hard spectrum of cosmic ray 
electrons a nd p ositrons as measured by Fermi [23, 111] 
and HESS [H H^l typically feature DM particles with 
masses on the order of c± 1 TeV or higher 13_, 31-33] . 

In many well motivated models, DM particles annihi- 
late with a cross section on the order of (av) ~ 10 
cm'^/s. If the dark matter's annihilation cross section (as 
evaluated in the early universe) is much larger than this 
value, the DM particles would have been overly depleted 
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in the early universe, resulting in a thermal relic abun- 
dance that is much smaller than the measured cosmolog- 
ical density. In contrast, in order for annihilating DM 
to account for the observed positron fraction, the DM 
particles must annihilate to leptonic final states with a 
cross section of {av) - 10^24 _ io-23 cmVs [slHIai]. 
Various proposals to accommodate this very high anni- 
hilation rate have been put forth, including annihilation 
cross sections which are enhanced at the low velocities 
found in the Galactic Halo (such as by Sommerfeld en- 
hancements) [l^, 111, I36l-l39] , DM particles which are pro- 
duced largely through non-thermal processes [4(| , or DM 
whose annihilation rate is highly boosted as a result of 
larger than expected inhomogeneities in their spatial dis- 
tribution. 

Further restricting any DM scenarios that might po- 
tentially account for the observed positrons is PAMELA'S 
measurement of the cosmic ray antiproton-to-proton ra- 
tio [41|, which is consistent with astrophysical expecta- 
tions. This observation strongly constrains the rate at 
which DM particles can annihilate to gauge bosons or 
quarks [III, |42, S^ ■ Furthermore, an annihilation rate to 
quarks, gauge bosons, or taus that is comparable to the 
leptonic rate required to produce the observed positron 
fraction would also lead to an unacceptably high fiux of 
prompt 7-rays from the Galactic Center |44l44a | and from 
the Inner Galaxy [43, 13 (constraints have also been de- 
rived from observations of dwarf spheroidal galaxies |49| - 
|52| | and the isotropic diffuse 7-ray background [53l45q ). 
As a result of these antiproton and 7-ray constraints, 
we are forced to consider DM candidates which annihi- 
late dominantly to electrons and muons, XX — ^ e^e^, 
XX — ^ /^^M I 01' to intermediate particles that later de- 
cay to combinations of e^e^, 11'^ ^j,^ , and tt+tt^. In con- 
sidering this latter case, known as eXciting Dark Matter 
(XDM) [Sg (see also Ref. |13|), we make use of the ana- 
lytic spectra described in Ref. [57| . 

High energy cosmic ray electrons and positrons un- 
dergo rapid energy losses as they diffuse through the in- 
terstellar medium. Above a few GeV, these energy losses 
are dominated by a combination of synchrotron and in- 
verse Compton scattering. To model the propagation of 
electrons, positrons, and other cosmic rays, we use the 
GALPROP package v54 (see Refs. [5l,|50|, and references 
therein), which includes up-to-date information pertain- 
ing to the local interstellar radiation field and the distri- 
bution of gas in the Galaxy. The former is relevant for 
the calculation of e^ energy losses via inverse Compton 
scattering, while the latter is relevant for the calculation 
of secondary cosmic ray species, produced in inelastic col- 
lisions of cosmic ray primaries with the interstellar gas. 
For synchrotron energy losses, we adopt a value of 5 /iC 
for the local magnetic field. By fitting both stable and 
unstable secondary-to-primary ratios in the cosmic ray 
spectrum, we can constrain the distributions of the inter- 
stellar gas and cosmic ray sources, as well as the timescale 
that cosmic rays reside within the Galaxy |60l l6l| . Codes 
such as GALPROP and DRAGON [62| assume a sim- 



ple diffusion zone with free escape boundary conditions 
(cosmic rays diffuse within the diffusion zone, but escape 
upon reaching any boundary of the zone) . We take this 
zone to be a cylinder, extending a distance L above and 
below the Galactic Plane, and radially 20 kpc from the 
Galactic Center. In this study, we vary the half-width of 
the diffusion zone over a range of values between L = 1-8 
kpc, in accordance with existing uncertainties [6l|. With 
further information from AMS (especially the measure- 
ment of the ^"Be/^Be ratio), this and other parameters 
of our propagation model will likely become significantly 
more constrained \63^. 

In Figs. [Hand [21 for the choice of L = 4 kpc, we show 
the cosmic ray positron fraction ($£+/($£+ -I- ^e-)) and 
the total lepton spectrum {E^dNe±/dE^±) predicted in 
a number of annihilating DM models. For each annihi- 
lation channel, we consider two values for the DM mass. 
The lower mass in each case is set to approximately the 
minimum value allowed by the new data from AMS (in 
particular, by the lack of any strong spectral cut-off). 
The upper mass in each case was chosen with the cos- 
mic ray electron spectrum (as measured by Fermi and 
HESS) in mind. As was found in previous fits to the 
PAMELA data, we find that annihilation cross sections 
in the range of ~ 4 x 10"^^ — 4 x 10"^'^ cm^/s are re- 
quired to accommodate the observed positron fraction. 
While changing the diffusion zone thickness or the radi- 
ation/magnetic field model can reduce this requirement 
to some extent, very large annihilation cross sections are 
a generic requirement of any annihilating dark matter 
model capable of generating the observed positron frac- 
tion. 

By considering different values of the diffusion zone 
thickness, L, we can slightly alter the predicted shape 
of the positron fraction, potentially enabling some DM 
models to better accommodate the measurements of 
AMS. In Fig. [31 we show for a representative DM model 
how changing this parameter impacts the shape of the 
positron fraction. Note that if we increase L (to 8 kpc, for 
example), we can soften the slope of the rising positron 
fraction, which in some cases allows for a better fit to the 
data. In Figs. U] and [5l we show predictions for various 
DM scenarios using choices of L = 2 and 8 kpc, respec- 
tively. As expected, the fits are worsened in the 2 kpc 
case, but are in some cases improved for L — 8 kpc. 

In none of the cases we have shown so far have we been 
able to find good agreement with the cosmic ray elec- 
tron-|-positron spectrum as measured by Fermi. To ac- 
commodate this measurement, we have to consider a cos- 
mic ray electron spectrum which is not a simple power- 
law, but instead breaks to a harder spectral index at 
high energies. This is not surprising, as at energies above 
^100 GeV energy losses limit the number of sources that 
contribute to the electron cosmic ray spectrum; the cos- 
mic ray electron spectrum at energies above a few hun- 
dred GeV is, therefore, likely to be dominated by only 
a small number of nearby sources. As a result, stochas- 
tic variations in the distribution of supernova remnants 
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FIG. 1: The predicted cosmic ray positron fraction (left) and electron+positron spectrum (right) in dark matter models 
annihilating to e^e~ (top) and to n'^fJ,' (bottom). The error bars shown represent the positron fraction as measured by AMS 
(black, left) and PAMELA (red, left), and the electron+positron spectrum as measured by Fermi and AMS-01 (black, right). 
In each case, we have adopted a propagation model that provides a good fit to the various secondary-to-primary ratios as 
described in the text, and with a diffusion zone half-width of L = 4 kpc. The expected backgrounds are shown as black dotted 
lines. For each annihilation channel, we show results for two masses. For annihilations to e+e~ and a mass of 350 GeV (900 
GeV), we have used a thermally averaged annihilation cross section of {av) = 4.4 x 10"^^ cm'^/s (2.5 x 10"'^* era? /s). For 
annihilations to /i^/i~ and a mass of 600 GeV (1.5 TeV), we have used a thermally averaged annihilation cross section of 
1.7 X 10"^'' cmVs (9.3 x 10"^* cm^/s). 



are expected to lead to local departures from the average 
cosmic ray spectrum found throughout the Milky Way 
(which may very well be a simple power-law) 26l |64| . 
When the Fermi electron-|-positron spectrum is taken in 
combination with the positron fraction as measured by 
AMS^ it is clear that more very high energy cosmic ray 
electrons are required than would be predicted using a 
simple power-law extrapolated from the low-energy spec- 
trum. 



In Fig. [SJ we show results using a broken power-law 
for the spectrum of electrons injected from cosmic ray 
sources, for the three DM models most capable of ac- 
commodating the observed positron fraction. Between 
4 and 100 GeV, we take the injected electron spectrum 
to be dNe/dEe oc E~^-^^ which provides a good fit to 
the observed low energy spectrum. Above 100 GeV, we 



harden this slope from -2.65 to -2.3."'^ With this spectral 
break, these three DM models now each appear to be 
capable of self-consistently accounting for both the mea- 
sured positron fraction and the overall leptonic cosmic 
ray spectrum. 



^ Note that instead of a hardening of the injected electron spec- 
trum above ^100 GeV, we could have instead considered a mild 
overdensity of local sources. As electrons from local sources ex- 
perience less energy loss than those from more distant sources, 
the contribution from local sources will exhibit a harder spec- 
trum, which could dominate the observed electron spectrum at 
high energies even if the shape of the electron spectrum injected 
from these sources is the same as that from the average source. 
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FIG. 2: As in Fig. 1, but for dark matter which annihilates into a pair of intermediate states 



which proceed to decay to 
e ' e (first row), to fJ,'^ fJ.~ (second row), to tv'^'jt~ (third row), and to a 1:1:2 ratio of e'^e~ , fi'^ fj,~ , and tt^tt" (fourth row). 
For annihilations to 2e"'"2e~ and a mass of 400 GeV (1.2 TeV), we have used a thermally averaged annihilation cross section of 
{av) = 7.3 X lO"^'^ cm^/s (6.2 x 10"^'' cm^/s). For annihilations to 2^+2^" and a mass of 800 GeV (2.5 TeV), we have used 
a thermally averaged annihilation cross section of 3.4 x 10"^'* cm'^/s (2.6 x 10~^^ cm^/s). For annihilations to 27r"'"27r~ and 
a mass of 1.0 TeV (3.0 TeV), we have used a thermally averaged annihilation cross section of 5.7 x lO"'^* cm^/s (4.1 x lO"'^^ 
cm^/s). And for annihilations to a 1:1:2 ratio of e'^e~ , /i'^ fi~ , and tt'^tt' final states with a mass of 500 GeV (1.6 TeV), we 
have used a thermally averaged annihilation cross section of 1.5 x lO"^** cm'^/s (1.3 x 10~^^ cm'^/s). 
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FIG. 3: The impact of the diffusion zone half-width, L, on the positron fraction in a representative DM model. In each case, 
we have chosen the diffusion coefficient to fit the non-leptonic background cosmic ray measurements. The dotted, dashed, 
dot-dashed, and solid lines correspond to L =1, 2, 4 and 8 kpc, respectively. The model used in this figure consists of a I TeV 
dark matter particle that annihilates to a pair of intermediate states which each decay to 7r"'"7r~, with a cross section given by 
16, 9.8, 5.7 and 3.5xl0"^^cm^/s, for L =1, 2, 4 and 8 kpc, respectively. 



III. PULSARS 

Pulsars are rapidly spinning neutron stars which 
steadily convert their rotational kinetic energy into radio 
emission, 7-rays, and cosmic rays, likely including ener- 
getic electron-positron pairs. When initially formed, typ- 
ical pulsars exhibit rotational periods on the order of tens 
or hundreds of milliseconds, and magnetic field strengths 
of ~10^^-10^'^ G. As a result of magnetic-dipole braking, 
the period of the pulsar's rotation slows down at a rate 
given by P = 3.3 x IQ-i^ (B/lQi^ G)^ (P/0.3s)-\ corre- 
sponding to an energy loss rate oiE = Att'^IP/P^ = 4.8 x 
lO^^erg/s (5/10^2 Q)2 ^p/Q 3 g)-4 (//io45g cm^). Young 

pulsars spin-down quite rapidly, typically losing a major- 
ity of their rotational kinetic energy in on the order of 
only ~10^ years. In both polar gap and outer gap models, 
a significant fraction of this energy can go into the pro- 
duction and acceleration of electron-positron pairs |65l - 

US- 

The spectral shape of the electrons and positrons in- 
jected from pulsars is often parametrized as |68| : 



^^E-"exp{-Ee/E,). 



(1) 



Although there is considerable uncertainty associated 
with these spectral parameters, a — 1.5 — 2.0 and Ec — 
80—1000 GeV cover the range typically found throughout 
the literature. We begin by calculating the contribution 
from the sum of all pulsars distributed throughout the 
Milky Way. To do this, we adopt the spatial distribution 
of the pulsar birth rate as described in Ref . [69| , normal- 
ized to an overall rate of 1 pulsar per century throughout 
the Galaxy, each with an average total energy of 10^^ erg. 
Following the procedure followed in the previous sec- 
tion, in Figs. [7] and [5] we show results for this distribution 



of pulsars. In Fig. [71 we adopt a simple power-law of 
index -2.65 for the injected spectrum from (non-pulsar) 
cosmic ray sources, while in Fig. [5] we break this spec- 
trum to -2.3 above 100 GeV. In each figure, we show 
results for two different choices of the pulsar spectral in- 
dex, a, and the diffusion zone half-width, L. In each case, 
we fix Ec — 600 GeV and normalize the positron and 
electron contribution from pulsars by assuming that 16% 
of the pulsars' total energy goes into high energy pairs. 
From these results (especially those shown in Fig. |S]), we 
conclude that for very reasonable choices of parameters, 
pulsars can provide a viable explanation for the observed 
cosmic ray positron fraction. 

In addition to the integrated contribution from all pul- 
sars throughout the Milky Way, there are two young 
and nearby pulsars which could each individually con- 
tribute significantly to the cosmic ray positron spectrum. 
The Geminga pulsar is 370,000 years old, 157 parsecs 
from the solar system, and pulsates with a period of 230 
ms. The pulsar B0656-I-14 (possibly associated with the 
Monogem supernova remnant) is considerably younger 
(110,000 years), although somewhat more distant (290 
parsecs), and more slowly rotating (P = 390 ms). These 
parameters, combined with their measurements of P, im- 
ply that Geminga and B0656-I-14 have each lost approx- 
imately 3 X 10^^ erg and 1 x 10^^ erg of rotational en- 
ergy since their births, respectively. If 4-5% of this en- 
ergy went into the production and acceleration of ener- 
getic e'^e" pairs, then these pulsars could be responsi- 
ble for the observed rise in the cosmic ray positron frac- 
tion [22, I23 . If we combine these two sources with the 
somewhat smaller contribution expected from the sum 
of all more distant pulsars [22], we estimate that if 3- 
4% of the total energy from pulsars goes into energetic 
pairs, this would be sufficient to account for the observed 
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FIG. 4: The same as in Figs. 1 and 2, but for a diffusion zone half-width of L = 2 kpc. For annihilations to 2e^2e~ and a mass 
of 400 GeV (1200 GeV), we have used a thermally averaged annihilation cross section of {av) — 1.2 x 10~ ^ cm /s (1.2 x 10~ 
cm^/s). For annihilations to 2^'^2fj,~ and a mass of 800 GeV (2.5 TeV), we have used a thermally averaged annihilation cross 
section of 5.6 x lO"^"* cm^/s (5.1 x 10"^^ cm^/s). For annihilations to 27r+27r" and a mass of 1.0 TeV (3.0 TeV), we have used 
a thermally averaged annihilation cross section of 9.8 x 10"^'* cm''/s (7.8 x 10~^* cm^/s). For annihilations to a 1:1:2 ratio of 
e^e~ , fi'^ fi~ , and tt'^tt' final states with a mass of 500 GeV (1.6 TeV), we have used a thermally averaged annihilation cross 
section of 2.7 x 10"^" cm^s (2.3 x 10"^^ cm^s). 
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FIG. 5: The same as in Figs. 1, 2 and 4, but for a diffusion zone half-width of L = 8 kpc. For annihilations to 2/i^2/i~ and 
a mass of 1.0 TeV (2.5 TeV), we have used a thermally averaged annihilation cross section of 2.9 x lO"'^* cm^/s (1.5 x 10~^^ 
cm /s). For annihilations to 2-k^2-k~ and a mass of 1.0 TeV (3.0 TeV), we have used a thermally averaged annihilation cross 



section of 3.5 x 10 cm'^/s (2.3 x 10 cm'^/s). For annihilations to a 1:1:2 ratio of e+e 



^+/i 



and TT^TT final states with 



a mass of 700 GeV (1.6 TeV), we have used a thermally averaged annihilation cross section of 1.6 x 10 cm /s (6.5 x 10 
cm^/s). 



positrons. 



IV. SUMMARY AND DISCUSSION 

In this paper, we have revisited both annihilating dark 
matter and pulsars as possible sources of the rising cos- 



mic ray positron fraction. Using the newly published, 
high precision data from AMS^ we have considered a wide 
range of dark matter models and cosmic ray propagation 
models. We find that models in which the dark mat- 



ter annihilates directly to leptons (e'' 



or [i^ [i 



are 



no longer capable of producing the observed rise in the 
positron fraction. Models in which the dark matter an- 
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FIG. 6: The same as in Figs. 1, 2, 4 and 5 but for a diffusion zone half-width of L = 8 kpc, and for broken power-law spectrum 
of electrons injected from cosmic ray sources {dN^- /dE^- oc E~^'^^ below 100 GeV and dN^^ /dE^- oc E~^'^ above 100 
GeV). The cross sections are the same as given in the caption of Fig. 5. With this cosmic ray background, the dark matter 
models shown can simultaneously accommodate the measurements of the cosmic ray positron fraction and the overall leptonic 
spectrum. 
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FIG. 7: The predicted cosmic ray positron fraction (left) and electron-l-positron spectrum (right) from the sum of all pulsars 
throughout the Milky Way, for an injected spectrum of dN^±/dE^± oc i5~l^'^^ exp(— J5g±/600GeV) and a diffusion zone half- 
width of L = 4 kpc (top), and for an injected spectrum of dN^± /dE^± oc iJ^^'®^ exp(— i5g±/600GeV) and a diffusion zone 
half- width of L = 8 kpc (bottom). For normalization, we have assumed that 16% of the pulsars' total energy goes into high 
energy electron-positron pairs. The error bars shown represent the positron fraction as measured by AMS (black, left) and 
PAMELA (red, left), and the electron-|-positron spectrum as measured by Fermi and AMS-01 (black, right). In each case, we 
have adopted a propagation model that provides a good fit to the various secondary-to-primary ratios as described in the text. 
The expected backgrounds are shown as black dotted lines. 
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FIG. 8: The same as in Fig.[71 but for a broken power- law spectrum of electrons injected from cosmic ray sources {dN^- /dE^- oc 
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nihilates into light intermediate states which then decay 
into combinations of muons and charged pions, however, 
can accommodate the new data (see Fig. [5]). In those 
dark matter models still capable of generating the ob- 
served positron excess, the dark matter's mass and anni- 
hilation cross section fall in the range of ~ 1.5-3 TeV and 



(av) - (6 - 23) X 10" 



-24 



Vs. 



We have also considered pulsars as a possible source 
of the observed positrons. In particular, we find that for 
reasonable choices of spectral parameters and spatial dis- 
tributions, the sum of all pulsars in the Milky Way could 
account for the observed positrons (see Fig. [8]) if, on av- 
erage, 10-20% of their total energy goes into the produc- 
tion and acceleration of electron-positron pairs (assuming 
a birth rate of one per century throughout the Galaxy, 
each with an average total energy of 10^^). It may also be 
the case that a small number of nearby and young pulsars 
(most notably Geminga and B0656-I-14) could dominate 
the local cosmic ray positron flux at energies above sev- 
eral tens of GeV. Taking into account these two excep- 
tional sources, we estimate that if 3-4% of the total en- 
ergy from pulsars goes into energetic pairs, these objects 
could be responsible for the observed positron fraction. 



Currently, we cannot yet discriminate between dark 
matter and pulsars as the source of the observed positron 
excess. We are hopeful, however, that future data fronr 
AMS may change this situation. In addition to contin- 
uing to improve the precision of their measurement of 
the positron fraction and extending this measurement to 
higher energies, AMS will also measure with unprece- 
dented precision a number of secondary-to-primary ratios 
of cosmic ray nuclei species, which can be used to con- 
strain many aspects of the underlying cosmic rays propa- 
gation model. Of particular importance is the ^"^Be/^Be 
ratio, for which existing measurements are limited to en- 
ergies below 2 GeV (kinetic energy per nucleon) , and with 
large errors (for a compilation of such measurements, see 
Tables I and II of Ref. (U). In contrast, AMS is ex- 
pected to measure this ratio with much greater precision, 
and up to energies of ~10 GeV. This information will 
enable us to break the longstanding degeneracy between 
the diffusion coefficient and the boundary conditions of 
the diffusion zone [63| . If these measurements ultimately 
favor propagation models with a somewhat narrow dif- 
fusion zone (L < 4 kpc), it would be very difficult to 
explain the observed positron fraction with any of the 
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dark matter models we have considered in this paper, 
while pulsar models could still provide a viable source 
for the positrons. Thus, in at least some plausible sce- 
narios, future data from A MS could allow us to rule out 
competing hypotheses and thereby finally reach a con- 
clusion as to the origin of the surprisingly large fiux of 
energetic cosmic ray positrons. 
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